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Abstract 

This technical paper of mathematical physics arose as an aftermath of 
Cassini's 2002 experiment [6], in which the PPN parameter 7 was mea- 
sured with an accuracy cr-y — 2.3 x 10^^ and found consistent with the 
prediction 7 = 1 of general relativity. The Orbit Determination Pro- 
gram (ODP) of NASA's Jet Propulsion Laboratory, which was used in 
the data analysis, is based on an expression (|8]) for the gravitational de- 
lay At which differs from the standard formula ([2]); this difference is of 
second order in powers of m - the gravitational radius of the Sun - but 
in Cassini's case it was much larger than the expected order of magni- 
tude m? /b, where h is the distance of closest approach of the ray. Since 
the ODP does not take into account any other second-order terms, it is 
necessary, also in view of future more accurate experiments, to revisit the 
whole problem, to systematically evaluate higher order corrections and 
to determine which terms, and why, are larger than the expected value. 
We note that light propagation in a static spacetime is equivalent to a 
problem in ordinary geometrical optics; Fermat's action functional at its 
minimum is just the light-time between the two end points A and B. A 
new and powerful formulation is thus obtained. This method is closely 
connected with the much more general approach of , which is based on 
Synge's world function. Asymptotic power series are necessary to provide 
a safe and automatic way of selecting which terms to keep at each order. 
Higher order approximations to the required quantities, in particular the 
delay and the deflection, are easily obtained. We also show that in a close 
superior conjunction, when b is much smaller than the distances of A and 
B from the Sun, of order R, say, the second-order correction has an en- 
hanced part of order iri^R/b^, which corresponds just to the second-order 
terms introduced in the ODP. Gravitational deflection of the image of a 
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far away source when observed from a finite distance from the mass is 
obtained up to 0{m^). 

1 Introduction 

In the framework of metric theories of gravity and the PPN formalism, 
the main violations of general relativity - those linear in the masses - 
are described by a single dimensionless parameter 7. The question, at 
what level and how general relativity is violated, in particular how much 7 
differs from unity, Einstein's value, is still moot. No definite and consistent 
prediction about it are available, except for the inequality 7 < 1, which 
must be fulfilled in a scalar-tensor theory, in particular those arising as the 
low-energy limit of certain string theories. To date, the best measurement 
of 7 has been obtained with Cassini's experiment, which has provided the 
fit (at 1-cr) 

7 - 1 = (2.1 ± 2.3) X 10"^ (1) 

Einstein's prediction is still acceptable, but more accurate experiment are 
needed and planned. 

While 7 controls also other relativistic effects, in particular those re- 
lated to gravito-magnetism, it mainly affects electromagnetic propagation. 
The differential displacement of the stellar images near the Sun historically 
was the first experimental effect to be investigated and is now of great im- 
portance in accurate astrometry. The bending of a light ray also increases 
the light-time between two points, an important effect usually named after 
its discoverer I. I. Shapiro fTP. Several experiments to measure this de- 
lay have been successfully carried out, using wide-band microwave signals 
passing near the Sun and transponded back, either passively by planets, 
or actively, by space probes (see [31], [24]). 

Cassini's 2002 experiment has implemented a third way to measure 

7 [3], in which coherent microwave trains sent from the ground station 
to the spacecraft (at that time about 7 AU far away) were transponded 
back continuously. The use of high-frequency carriers (in Ka band, 34 
and 32 GHz) and the combination with standard X-band carriers (about 

8 GHz) allowed successful elimination of the main hindrance, dispersive 
effects due to the solar corona traversed by the beam. The tracking was 
carried out around the 2002 superior conjunction; the minimum value 
of the impact parameter of the beam was 1.6 ii0, but in effect only 18 
passages have been used, with a minimum impact parameter of ~ 6 Rq. 
The two-way total amount of phase between the time of emission and the 
time of arrival has been continuously measured in each passage. In effect, 
however, NASA's Deep Space Network provides the phase count in a given 
integration time r. Mathematically, in the limit r — » this would give the 
received frequency, in which Doppler effects and gravitational frequency 
shift are mixed up (Sec. 4). Cassini's observable, therefore, can also 
be assessed in terms of the predicted change in frequency, as in [3]; but 
in practice, taking r small would introduce unacceptable high-frequency 
noise. The change in light-time in a given integration time is the correct, 
theoretically available observable. 



2 



In the standard formulation for a superior conjunction, and taking the 
Sun at rest, the (one-way) light-time from an event A to an event B is: 



ts-tA= r,, +At^ r,, + (1 + 7)mln 1 1 ■ (2) 

'a ' 'b ^ab 

where m = 1.43 km is the gravitational radius of the Sun, rj^,rg are, 
in Euclidian geometry (See Fig. 1 left), the distances of A and B from 
the Sun and r^^ their distance. The velocity of light c is unity. At, the 
increase of the light-time over r^^ , is the gravitational delay. 

In a close superior conjunction A and B are on the opposite sides of 
the mass and the Euclidian distance 6o of the straight line AB from the 
mass fulfils bo <^ {'Ta^'^'b) ~ 0{R), say. In this approximation eq. ([2]) 
reduces to 

tB-iA=r^^+At = r^^+{l+j)mlnl^-^^j , (3) 

with a logarithmic enhancement over the formal order of magnitude At = 
Taking the logarithm equal to 10, this provides an estimate of 
the timing accuracy in terms of the error in 7: 

o-^j = 1.43 o-^ X lO'^cm, (4) 

corresponding, in Cassini's case, to 30 cm. (|3]) embodies also the one-way 
frequency change Av induced by gravity between A and B. Their motion 
makes 60 (and the distances) change with time, so that, for a one-way 
experiment. 

The basic geometric setup is straightforward: a point mass m at rest 
at the origin in an asymptotically flat space generates a line element with 
rotational symmetry. An invariant Killing time t is defined; events on 
each t = constant surface are 'simultaneous' and the metric components 
are constant. The proper time ds — \/ goo (r) dt of a static observer differs 
from dt by the red-shift factor y^gooir)- A null geodesic runs from the 
event A (with radial coordinate r^ and time t^ ) to the event B (with radial 
coordinate and time t^); it stays on a plane, taken here as the equato- 
rial plane 6 = -k 12. The (invariant) longitude difference = cj)^ — (j>j^ 
completes the setup. In the PPN formalism and isotropic coordinates the 
metric reads: 



da^ = A{r)dt^ - B(r)df = 




where 

^As stated in the supplementary material, in eq. (2) of [B] the two terms in the right-hand 
side should obviously be multiplied by a factor 2. This error, of course, had no consequence 
on the computer fit. 
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= dr^ + r'^{de^ + s\aed(t>^) = dr^ + r^dfi^ 



is the Euclidian line element. The parameters 7, /3 and e are equal to 1 
in general relativity; while 7 and (5 are accurately known, currently no 
information is available about e. 

In our case the best mathematical tool to deal with electromagnetic 
propagation is not null geodesies, but the theory of eikonal. It is known 
(e.g., [T^) that in this problem Fermat's Principle holds, corresponding 
to the refractive index 



we develop ab initio the eikonal and solve for it by separation of variables 
(Sec. 4). The radial part provides Fermat's action as a radial integral 
containing N{r) and the impact parameter h; when computed at the true 
value /itruo, such action is just the required light-time. The solution can 
be obtained recursively, using appropriate expansions in powers of m: 
the expansion for h begins with ho = bo, the distance of the straight 
line AB from the origin. In this way the variational nature of the problem 
brings about a great conceptual and algebraic simplification. At the linear 
approximation in m one would expect that the light-time contains h\ , the 
correction in the impact parameter linear in the mass; as one can see from 
(12]), this is not the case. This property is generally true: the correction to 
the light-time 0(m)* does not contain /i^ (Sec. 6). 

Cassini's and many other space experiments have been analyzed us- 
ing NASA's Orbit Determination Program (ODP), developed by NASA 
at Jet Propulsion Laboratory in the 60's and steadily improved since; a 
new version called MONTE is under development. The ODP, whose the- 
oretical formulation is due to T. D. Moyer [2T], integrates the equations 
of motion of the relevant bodies and provides their trajectories in the 
ephemeris time. This task is carried out in a reference system - called 
BCRS (Barycentric Coordinate Reference System) - in which the centre 
of gravity of the solar system is at rest and the Sun moves around with 
a velocity vq ~ lOm/sec = 3 x 10~*c. As discussed in [2], the light-time 
in this this frame difi'ers from the rest frame of the Sun essentially due to 
Lorentz time dilatation; being of order vq, this difference is quite below 
the sensitivity of Cassini's experiment. We do not discuss this point any 
more; t is just Killing time. 

The ODP uses a fictitious Euclidian space Sz{x,y,z), which corre- 
sponds to the isotropic coordinates of (|6}. This space is just a computa- 
tional convenience and should not be considered as a physical background 
in which gravity acts. For example, replacing r, the Euclidian distance 
from the origin, with r + km, where k is an arbitrary constant, is fully 
legitimate in a covariant theory, but it destroys the conformal fiatness of 
space, introduces a gravitational potential —krr^jr^ and adds a second- 
order term to the delay At. Strictly speaking, the word 'delay' is inap- 
propriate: we just have a light-time and there is nothing with respect to 
which a delay can be reckoned. The object of the measurement is the time 
change of the delay. The arbitrariness of the radial coordinate also affects 
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gravitational bending: its second-order approximation up to 0(m/}if , 
depends on which radial coordinate is used (see [TH], [10] and [55]) [7]). 

It should also be noted that the spacetime coordinates of the end 
events are not directly provided in the experimental setup and depend on 
the gravitational delay At, the very quantity one sets out to measure. The 
trajectories {t) and (t) are given by the numerical code; the starting 
time tj^ is just a label of the ray, but the arrival time is greater than 
tj^ + r^g . The way out is to take for the end point 

where = dr^/dt For a typical velocity 10~* c the correction is of order 
20 X 1.4 X 10^ X 10"* = 300 cm, and the a priori accuracy in At is 
sufficient. 

Since for electromagnetic propagation dt and d£ in ([6]) are almost equal, 
([2]) is the correct approximation to the delay to 0(m); one would expect 
this to be the first term in an expansion in powers of m/bo, so that the 
next term should be 

m- m Rq „ q ^© 
m— — m— ; — — 0.3 -; — cm, 

&o rt© 00 Oo 

quite below Cassini's sensitivity. The present paper arose because the 
ODP (eq. (8-54) of [21]), in fact does not use ([2]), but, in our notation. 



' AB 



+ At = + (1 +7)mln 



■^A +rAB + (1 +7)w^ 

^A - '•as + (1 +7)"i 



(8) 

We have not been able to fully reconstruct Moyer's derivation of this 
expression. It introduces non linear corrections arising from non linear 
effects of linear metric terms, but no quadratic metric terms. However, 
the difference between the two expressions of the delay is much larger 
than the estimate above; this arises because in Cassini's case, in Q the 
denominator rj^+r^— r^^ is much smaller than the numerator ~ 2r^g . 
Indeed, 



At - (At)„,, = -2(1 + 7)^^-^ = -(1 + 7)^^, (9) 
where we have introduced the harmonic mean of the distances 

1^^ + ^^^. (10) 

R r^r^ 

If, as in Cassini's experiment, r^ ^ rj^ = lAU = 200 ii©, R = 400 7?© 
the correction is about 

1600 m^(':?0V = 500 ^^°'' 



Rq \ bo J \ bo 

Even at ~ 6 Rq this correction is somewhat below the sensitivity (|4} and 
it should not have affected the result. However, it cannot be excluded 
that neglected non linear terms relevant for Cassini's experiment affect 
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the fit U]). One could say, ((8)l is mendacious; a full clarification of the 
problem is needed. 

Empirically dropping or keeping 'small' terms may lead to inconsisten- 
cies and does not work; the rigourous method of asymptotic perturbation 
theory (see, e. g., [TT], [TS]) must be used. We briefly sketch it now at a 
practical level. One begins with a wise choice of a dimensionless 'small- 
ness' parameter, and expands every function in the corresponding power 
series. Our main choice will be m/bo, but convenience may suggest using 
other lengths, like in m/r. An asymptotic series 

is a formal object assigned just by the sequence of its coefficients Gs; 
arithmetics and calculus follows the obvious rules for sum, multiplication 
and differentiation. Equality between two asymptotic series just means 
that the coefficients of the same order are equal. The value of G{m) as 
a function of m does not play any role, and even the convergence of the 
series is irrelevant; what matters is only the truncated value at any order 
k 

fc / \ s / \ fc + 1 

s — 

The parameter should not be understood as a fixed number, but as a vari- 
able which tends to zero. The symbol 0{.) means order of infinitesimal; 
it states how fast the remainder tends to zero as the parameter diminishes. 
An asymptotic series can be constructed from an ordinary arbitrary func- 
tion G(m); but a whole class of functions give rise to the same series; for 
example, if Ga is the sequence generated by G{m), the same sequence is 
also generated by 

[l + Pexp(-Q6o/m)jG(m) (g > 0). 

In this way any recursive iteration then proceeds automatically and safely, 
even in the most complex situations. 

In our case light-time will be provided as an asymptotic power series 

+ A. (ni, 

with dimensionless coefficients As . Ai provides the lowest, standard 
approximation to /S.t (see ([2|. In principle, asymptotic analysis does not 
provide a numerical estimate of the remainder in a given situation; this is 
a physical, not a mathematical question. But when the problem, properly 
formulated, does not contain small dimensionless quantities other than the 
smallness parameter itself, one can expect the mathematical operations 
leading to the result to maintain the order of magnitude and to lead to 
expansions whose coefficients are numerically of the same order. This is 
the case of deflection, the angle between the asymptotes of the ray. There 
is only one length in the problem, the distance b of the point of closest 
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approach, or, equivalently, the impact parameter h = bN{b) (see Fig. 4); 
hence in the expansion 

s 

the coefHcients Ss are dimensionless numbers, solely determined by the 
PPN parameters and, must be of order unity (see Sec. 9). But in the 
delay problem the coefficients As depend on the geometrical configura- 
tion. They are of order unity in the generic (but scarcely interesting) 
case in which r^,r^ and bo are of the same order; but in a close supe- 
rior conjunction - of crucial relevance in experimental gravitation - when 
bo ^ (i'aj^b) — 0{R), besides m/bo, there is another smallness parame- 
ter, namely, bo/R, and there is no reason to exclude that the As increase 
with R/bo beyond the expected order of magnitude unity. This we call 
enhancement. We already saw in (|3} that Ai is enhanced, albeit only 
logarithmically; the ODP correction formally of second order, is en- 
hanced by R/bo- This could place serious limitations on the method and 
even invalidate the iteration itself. This would occur, for instance, when 
mR « bo; if 60 = ^^0 = 1/200AU, this corresponds to = 2000 AU. The 
enhancement, which has never been discussed in the literature, has been 
fully understood and tamed in the present paper (Sec. 8). We have found, 
indeed, that the second-order terms embodied in the ODP expression 
which was used in Cassini 's experiment are just the enhanced second-order 
terms; Cassini's result is still safe. 

The problem can be reduced to one of ordinary optics; due to its varia- 
tional nature, the eikonal function can be easily solved in an expansion in 
powers of m/h. The second-order expression of the light-time for a static 
spacetime has been obtained; extension to third order is also easy. This 
approach should be compared with the much more general work of [18] , 
who consider Synge's world function n(a;^,a;^) in a generic spacetime 
for a generic geodesic (not necessarily null) between two events A and 
B. On the basis of Hamiltonian theory, they develop a method to solve 
for Q,{xj^,Xg) in a formal power series with respect to the gravitational 
constant G and compute it up to the second order. In the null case the 
world function vanishes on the solution and becomes the eikonal function. 
Out method, limited of course to the spherically symmetric case, exploits 
directly the variational nature of the problem and leads to the second- 
order expression of the light-time, which agrees with the expression of 
[18] : extension to third order is also easy. 

For a realistic observation of a distant source from a point B at a finite 
distance r^ , ()13[) must be generalized to an expansion of the type 

*.=E*..(^)(^)-. (") 

s 

where h is the impact parameter. The linear term has been evaluated in 
[19] , $ 40.3; the quadratic correction will be obtained in Sec. 9 . 
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2 Hyperbolic Newtonian dynamics 



Newtonian dynamics of a test particle attracted by a point mass AI, an 
exactly soluble problem, illustrates these issues. We consider a motion 
in the equatorial plane 9 = 7r/2, with radial coordinate r and azimuthal 
longitude 0. The Lagrangian function 



r _ 1 

A-Ncw — — 



dt \dt 



GM 
+ — 



r 

keeps the total energy constant; Voo, the ultimate speed of the par- 

ticle at a large distance, plays a role analogous to the speed of light and 
will be taken equal to unity. Then 

where m = GM/v'^ is the gravitational radius. is an ignorable coordi- 
nate, so that the angular momentum 

d{d(l)/dt) dt ^ ' 

is constant. Since the velocity at infinity is 1, ft is also the impact param- 
eter. Eliminating dt we get: 

r^ = ± , ^ ; (18) 

hence 

h = ^/b{b + 2m) (19) 

determines b, the distance of closest approach where dr/d(t> — 0. The sign 
depends upon whether the ray is ingoing or outgoing. Integrating we get 
the true anomaly 

/ = arccos ^f!"'^ J'^^ ^ ■ (20) 

Alternatively, the motion can be expressed in terms of the semi-major 
axis a = m and the hyperbolic eccentricity e = 1 + b/m: 



a{e^ - 1) 
r = — — . 

1 + e cos / 

The acute angle S between the asymptotes is given by 



(21) 



sin (5 = sin f 2arccos ( — - | | = t—— — i / 1 — — ttt (22) 

V Veyy h + m\ (b + m)^ ^ ^ 

This angle has a regular expansion in powers of m/b, with no enhance- 
ment. 

Consider, however, the hyperbola determined by two points A and B 
on the opposite sides of the vertex (right side in Fig. 1). As in space 
navigation - in particular in the ODP - the end points are provided in 
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Figure 1: The background Euclidian geometry. The mass at O and the end 
points at A and B define a triangle AOB; the distance OH = bo from the 
straight line AB to the mass at the origin is taken as the unit of length. 
The angles a are taken positive. The internal angle can be obtuse 

(left) or acute (right); in the first, more interesting case, when, in addition, 
— 1 ^ ^^0: wc have the most important case of a close superior conjunction, 
in which the deflection is large. Elementary trigonometry gives the relation 
6o +r2 - 2r^r^ cos = sin . 
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terms of the initial and final position vectors, or equivalcntly, in terms of 
the initial and final distances and and the elongation angle "5^^ ; 
the "unperturbed distance of closest approach" may then be calculated 
from elementary geometry: 

■ _ r ^r^ sin$^^ 

V^A +'"1 - 2r^'"B cos*^^ 
Choosing the angles , and positive, we can express b in terms 
of 6o with the condition 



■, bo bo 

^AB ~ "^A + '^B ~ arccos h arccos — = 

''"a ^ b 

/b(6 + 2m) m \ 

= arccos -r — + 

\r^(b + m) b + mj 

+ arccos f^(^±^ - -^V (24) 
\rg{b + m) b + mJ ^ ^ 

The symmetric case r^^ = — R is sufficient to exhibit the problem. 
The condition reads: 

R b + m b + m R ^ ' ^ ^ 



6^ + (2m- 1)6 - m(i? + 1) =0, (26) 



with the solution 



2b bo-m^ L^mRpn^ ^^7) 



bo bo V bo bo 

Expansion in powers of m gives 



Vo-'^[vo-')ro-[[b^J "VUJ • ^''^ 

The enhancement is clear: when R — O{bo) the truncation error at order 
k is 0(771/60)*^^^, with a coefficient of order unity, as naively expected; but 
when - as in a close superior conjunction - R ^ bo, the error is larger, 
0{mR/bl)''+\ Formally this requires introducing another smallness pa- 
rameter 60 /R and expanding every coefficient of the primary m-expansion 
in descending powers of R/bo- Of course, the condition 

^ « 1, (29) 
"0 

must be fulfilled, lest the whole procedure breaks down. One could say, 
anchoring the trajectory at far away end points has a lever effect, so that 
an increase in the mass produces a large increeise in closest approach. 

The quantity (|29p gives, in order of magnitude, the ratio between 
the deflection ~ m/bo and the angle bo/R which separates the central 
mass and a distant star, as seen from a distance R. Hence the limiting 
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constraint above implies that the geometry of astronomical deflection is 
the same as in the classical case (see Fig. 4): sources in the sky near the 
Sun are displaced outward by an amount inversely proportional to the 
angular distance. The transition through the milestone mR — marks 
the passage to the gravitational lensing regime, in which the image can 
appear on both sides. 

In Sec. 8 the light-time enhancement is dealt with in the general 
case and it is shown that the dimensionless coefficients As in (I12[) are 

o{R/bor-\ 



3 The radial gauge 

The metric of a spherical body at rest has the general form 

ds^ = A{r)dt^ - B{r)dr^ - C{r)r^dQ.'^, (30) 

where A{r), B(r),C(r) are of power series the form: 

m 
r 



It is asymptotically flat, so that Ao = _Bo = Co = 1. The radial coordinate 
is otherwise arbitrary; this is the gauge freedom at our disposal. For 
consistency, however, any change r ^ f = g{r) must become an identity 
at infinity and have a similar expansion: 

g{r)^r + gim + g2— + (32) 

the coefficients As,Bs,Cs axe not gauge invariant. Two gauges are com- 
mon. In the isotropic form - the canonical choice in space physics - 
C{r) = B{r), so that 



ds^ = A{r)dt^ - B(r){dr^ + r^dO.'^) = A{r)dt^ - B{r)df; (33) 
the space part of the metric is conformally flat. We define 



In the PPN scheme (e. g., [31]) 

In 'Schwarzschild' gauge C{f) — 1 and 

ds^ = A{f)dt'^ - B{f)df^ - f^dn^; 

the area of a sphere of radius f is just the Euclidian expression 47rf^, 
which deflnes f in an invariant way. In the original Schwarzschild solution 
A{f) = 1/B{f) = 1 — 2m"// f. To get the isotropic form one requires 
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g\r)^B{g{r))(^y r'; (36) 



to first order 



r = r + jm + . . . . (37) 
In the present paper a third radial coordinate 



p = rN{r)=r^^^ = r + mNi+m^^ + ... (38) 

plays an important role. It is a monotonic function of r and ensures 
A{p) = C{p). In the linear approximation it was introduced by Moyer 
in [21] (eq. (8-23)), and boils down to just adding to r a constant term, 
equal to 2.95 km for the Sun. 



4 Geometrical optics 



It is convenient to reduce the problem to geometrical optics using the 
eikonal function 6. In a generic spacetime 6 fulfils the eikonal equation 

g'"'^^,e^,& =0; (39) 

its characteristics are the null rays (see, e. g., T). © is the phase of the 
electromagnetic wave. Let = r''(s^),r^ — r'^{Sg) be the trajectories 
of the end points, given as functions of their proper times s^^jS^; let 

ds^' dSg 

be the corresponding four- velocities. Clocks associated with them measure 
the proper frequencies 

= -<a,6 = g, =-<a,.6 = ^. (40) 

In the simple case in which the end points are far away from the source, 
where the metric corrections can be neglected, the contribution to the 
frequency difference corresponds to the ordinary Doppler effect, and can 
be evaluated with a slow motion expansion; the change in S between A 
and B is determined by the accumulated gravitational effect along the ray 
and mainly come from the region near the mass. 

ff^'a^ea.e = o = iv^(r)(9t6)^ - v© ■ ve, (4i) 

where V is the Euclidian gradient operator. We are really interested only 
in the spherically symmetric case, but the reasoning of this Section holds 
also for an arbitrary N{r). 

& is the phase; propagation occurs keeping it constant. Separating 
space and time variables with 

6 = 6t(t) + 6r(r), 
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leads to the class of solutions 



G=u:o[&{r)-tj, (42) 

where cJo©(r) is the spatial part of the phase and ujq is a constant fre- 
quency. & has the dimension of time and satisfies 

Ve • V6 = iV^(r). (43) 



If a clock is at rest relative to the mass, v'^ — (1, 0)/ y'Ajr), and the mea- 
sured proper frequency wo / \/ A{r) includes the appropriate gravitational 
shift away from the asymptotic value ujq. This is enough to reduce the 
problem to geometrical optics (see, e. g., [^, Ch. III). A ray r(^), as func- 
tion of the Euclidian arc length £, is orthogonal to the eikonal surfaces 
©(r) — const and fulfils 

^N{r)^)=\7N{r). (44) 



d£ V d^. 

The index of refraction is the rate of increase of the spatial phase along 
the ray: 

Consider now Fermat's action functional 



SHX)] ^ £ ^AiV(r)/f^ ^ £ dXC^, (45) 

where the trajectory, any path joining the end points, is expressed in 
terms of a generic parameter A: 

r(AJ=r,, r(A^)=r^. (46) 

Since the action is, in fact, independent of the choice of A, no generality is 
lost if dX = d£, the Euclidean line element. The Euler-Lagrange equation 
for the action (|45ll reduces to (|44[l . The actual elapsed time 

t^-t^^ S{A, B) = dlN{r) = - 6^ (47) 

is just the value of S[] computed at a local minimum - the actual ray 
{Fermat's Principle). One should keep in mind the distinction between 
the action functional, with its argument in square brackets, and the ac- 
tion computed at the extremum, an ordinary function of the end points 
denoted with S{A,B). In S{A,B), but not in S[], it is allowed to replace 
the generic independent variable A with a more convenient one related 
to the solution, like r. For simplicity, the different functions denoted by 
the symbol S are distinguished by their arguments; below, the quantity 
5(r^ ,rg\b) = S{h) will be introduced to denote the action corresponding 
to a ray anchored at and r^, but with arbitrary b (or h). 
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5 The solution 



The eikonal function provides a deep simplification in the evaluation of the 
light-time. Having already separated out the time, the three-dimensional 
eikonal equation (|43|l in spherical symmetry and in the equatorial plane 
can be solved by separating out the longitude (j): setting 

e(r,0) = 6,(r)-h©4<^). 

It satisfies 



so that ©0 is a constant. Setting = hcf), the eikonal equation reduces 
to 



with the primitive 

^Vr2iV2(r)-ft2. 

The -I- and the — signs correspond, respectively, to an outgoing and an 
incoming photon. The radial coordinate of closest approach 6, where 
©r = 0, is the solution of 

hN{b) = h; (48) 

since r > 6, © is a real function. In Sec. 9 it will be shown that h, just 
like in the Newtonian case, is the impact parameter (Fig. 4). The total 
phase is, therefore, 

&=uo{h<j)±j ^^/r'^N^ir) -h^ -t^ . (49) 

A wavefront propagates keeping © constant, so that the time along the 
ray is 



t = ±l —^Jr'^m{r)-K^ + hct). (50) 



In the usual case (see Fig. 1), in which the angle AOB is obtuse, the 
ray has two branches, both taken with the positive sign: an incoming one 
from to b and an outgoing one from fo to . In the acute case b is 
never reached and we have just an outgoing ray from to . In both 
cases, in going from A to B the longitude increases by (p^ ~ (f)^ — <I>^^ . 
The quantity 



S(h) = f —Jr^mtr) - /i2 ± ^^r2Ar2(r) - + h^.„ (51) 
Jb r Jb r 



For a function of a single variable a prime indicates the derivative. 
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Light-time 




True h 



h 



Figure 2: The minimum of the reduced action ((5T|) is equal to the light-time at 
the true value /itruc- 

gives the phase change, hence the light-time, between the end points, but 
the quantity h is still arbitrary. The upper (lower) sign corresponds to 
the case in which the angle AOB is obtuse (acute); in the latter case the 
two integrals combine in a single one from to , and h disappears as 
a lower limit. (|5H) is what Fermat's action functional becomes when its 
variability is restricted to h and the longitude constraint is not imposed; 
it shall be called reduced action. At the true value it satisfies 



keeping the end points fixed. 

The present work aims at providing the theoretical foundation for the 
time delay in all configurations; the sign freedom allows dealing with both 
cases at the same time, but applications will be mainly given for a con- 
junction, with the The origin of longitudes is arbitrary. This general 
approach is relevant, for ex ample, for a spacecraft on an almost parabolic 
orbit, as in the Solar Probe concept; with a perihelion as low as 4_Rq, it 
can have a strong enhancement of the light-time even in the acute config- 
uration. 

In the derivative S'{h) there are no contributions from the lower limits; 
then (|52p provides h as an implicit function of the total total elongation 
$Ai3: 



5'(;itruc) =0, 



(52) 




(53) 
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Hence ((51)) readfl 



f-^B rN(r) 
S{h)^ drN{r) ± / driV(r 



rN{r) 



lb \/{rN{r)y ^ Jb \/{rN{r)y^ -K'' 

(54) 

i 

Both integrals are convergent (and in the acute case the singularity 
at rN{r) = h is not even reached). (|51|l suggests the introduction of the 
function 

/■r fir , 

G{r,h)^ —sJ{rN{r)Y~h^, (55) 

Jb ^ 

in terms of which 

S{h) = G{rB,h)±G(rA,h) + h^AB. (56) 

dsn readfl 

Gh (rs , ft) ± Gh (rA ,h) + ^AB^O. (57) 

While in (|51|l h is an independent parameter, in (|54p it is fixed by (|53p . 

This expression for ft can also be derived directly from Fermat's Prin- 
ciple, thus providing its significance. Fermat's action (|47|l . expressed as a 
function of r, has the Lagrange functional 



CA<f>{r)] = N{r) ^l + r\d<t>/dr)\ (58) 
with the (positive) constant of the motion 

=±^^a^^ = ft. (59) 



d{d(t>/dr) y/TTrHd^f/drf dr 
The upper (lower) holds for the outgoing (incoming) branch. Integrating 

d<t> , h , ^ ian\ 

r— = ± , = ± , , (60) 

dr ^r2iV2(r)-ft2 _ ft2 



(|53[l is recovered. Comparison with the Newtonian case (|18p shows that 
the latter corresponds to the exact index of refraction 



/ 777 

iVNcw(r) = Jl + 2--, (61) 

corresponding, as expected, to 7 = 0, A'^i = 1 and N2 — —Na — —1/2, 
etc. 



•^In a sliglitly inconsistent notation, we often use h to denote both an independent and 
variable quantity, and the fixed value fttruc determined by the elongation. The context should 
be sufficient to clear the ambiguity. 

*The suffix indicates partial derivative. 
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Figure 3: Three ways to define the separation of the ray from the origin: the 
distance bo = ho = OH (in this paper often taken as unit of length) of the 
straight hne AB; the distance b = OK the point of closest approach; the impact 
parameter h = bN{b) = OJ . 
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6 A variational argument 



At this point one could proceed as follows: using power series, solve (|53p 
for h in terms of $as, a known quantity. The value of h, inserted into 
(|54p. provides the required light-time. The stationary character of the 
action H52p . however, brings about a deep and important simplification. 
This is already tacitly applied in the usual derivation of the gravitational 
delay ((2)1. To first order, the integral of dt — N{r)d£ in p3[) reads 

tB-tA = / de+ij + i) / de-; 

the second integral can be carried out along the straight path from A to 
B, leading to the characteristic logarithmic term. In principle, however, 
the first integral should take into account the (first order) deflection; we 
should understand J dl = as the Euclidian length of the bent arc 
between A and B. But the length r^^ of the straight segment AB is a 
minimum in the set of all curves joining A and B, so that l^^ — r^^ 
vanishes to 0(m)lf| Ray bending is irrelevant here. 

In order to exploit the variational nature of the problem it is convenient 
to apply power expansions before imposing the extremum condition (I52|) . 
We just need the value of the reduced action (|5ip S{h) — "^Zs^^^^i^) 
at the value which fulfils ((52}, namely, = m" S's{h). Setting h = 
ho + mhi + m^/i2 and expanding, the solution to second order is obtained 
iteratively: 



S'oiho) = 0, (62) 
/iiS;'(/io) + S[iho)=Q, (63) 

/i2S^,'(/io) + -^S^;'{ho) + hiSUho)+S'2{ho)^0. (64) 



In the expression 



Sih) = So{ho) + m{hiS'oiho) + Siiho)) + 

+ (^h2S'o{ho) + ^^S^{ho) + hiS[{ho) + S2(ho)J , (65) 

the effect of the extremum property is clear: since S'o{ho) = 0, the first or- 
der term does not contain hi, and the second order term does not contain 
/12; in general, the term in S{h) of order m'' does not depend on h^. This 
important result is refiected in the general approach of [18]. Referring 
to the equation numbering of that paper, their world function i}{Xj^,Xg) 
fulfills the Hamilton- Jacobi equation (30). In the null case SI = 0, (30) 
becomes the eikonal equation. Their Theorem 2 proves that the ri'^- 
order Q," can be expressed in terms of integrals along the lowest order 

didactical remark is in order here. This minimum property, crucial to the argument, is 
often omitted in the usual derivation. See, e. g., 1191 p. 1107, |9] p. 125; in equation (17.59) 
of [5], p. 581 the minimum is not mentioned and a factor 4 is missing in the argument of the 
logarithm. 
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Minkowskian path. In our case this variational Lemma clarifies the mat- 
ter and produces considerable simplifications. Using (|63|l . the light-time 
to second order reads: 



S(h) = So{ho)+mSi{ho)+m'' i^-jSo{ho) + hiS[(ho) + S2iho)j , (66) 
where hi is given by (|63p . The delay coefficients (|12|l read 

Ai = Si(fto), A2 = Y^o(^o) + hiSiiho) + S2{ho). (67) 

The second-order correction in the impact parameter ^2, given by (|64p . 
is needed only at third and higher orders. For the record, note the third- 
order contribution to the light-time: 

A3 = hihiS'o + -^S'o + h2S[ + -j-S'l + hiS2 + S3, (68) 

where, for simplicity, the arguments ho have been understood, and /12 is 
provided by (|64p . 



7 Power series 

We now proceed to apply this simple and general Lemma to the light-time. 
To lowest order, in we use b = bo — ho and N{r) = 1, so that 



Soiho) = ^rl~hl±^rl-hl + 

— ho ( arccos — ± arccos — — <I> ) . (69) 



The condition S'o{ho) — determines ho with the trigonometric relation 
(see Fig. 1) 

$ = arccos — ± arccos — , (70) 

' R 'A 



''^'■«.sin$^^. (71) 



Therefore 



So{ho) = -hl± ^r\-hl = r = ^ + rl - 2r ^ cos$^^, 

(72) 

is the geometric distance AB. ho is now fixed and can be taken equal to 
unity without loss of generality. Because of the variational Lemma, at the 
next order we can retain h = ho = 1; (fSTjl . with N{r) = 1 -|- mNi/r, reads 
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In the obtuse case (with the + sign) the logarithm has the argument 

as easily checked by cross multiplication using (|72p : then the standard 
expression ((2]| is properly recovered. (See the Appendix for the confusion 
that can arise due to the gauge freedom and the difference between closest 
approach and the distance bo). In the acute case, instead, 

Before proceeding to the next order we need to evaluate hi with (|63p . 
Differentiating (f69)) twice we easily get 



1 



1 



(75) 



In Sec. 9 the obtuse case in which oo will be considered; it simply 

gives 

hi=Nir^. (76) 
Considerable simplification may be achieved with the aid of the identities: 



'% -hl = rg{r^ 



cos'1>^b)Aab; 



±rjr^ - Tg cos$^^)/r^ 



In both cases the expression for h\ becomes 



fti =iVi 



1 — cos <E>^ 
sin $ ^ „ 



(77) 
(78) 

(79) 



It is useful to record the value of h\ — hi — Ni : 



Enhancement is at work: in the obtuse case, with the + sign, the elonga- 
tion comes close to tt and hi becomes large, as discussed in the following 
Section. In the acute case hi remains of order unity. 

At the next order (see (|65p i. we need G2{r,h), Gi{r,h) and its first 
derivative with respect to h, and Go{r, h) with its first and second deriva- 
tives (see (|55p ). At order s we need Go{r,h) with its first s derivatives. 



bi 
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If these difFerentiations are carried out before the integration, a technical 
difficulty arises, h appears both in the lower limit and in the square root. 
In the obtuse case, already at the second order each of the two contri- 
butions diverges; the second derivative of the integrand, for instance, has 
a non-integrable term cx (r^ — /iq)"'^''^; it turns out, however, that this 
divergence is compensated by the lower limit contribution. At higher or- 
ders the complexity increases. In the acute case the singular point is not 
within the integration domain and no hindrance arises. This suggests that 
the integration is best carried out first, leading to a finite result whose 
differentiation is straightforward. 

The hindrance arises because as m ^ the singular point at r = fe 
moves. The integration variable 

keeps the singularity fixed ai u = 1 and cures the problem. Note the 
appearance of Moyer's radial coordinate 

p = rN(r) = r + mNi+7n'^ — . (82) 

r 



Then ((55| reads 



"("■^ dlnr(w) 



G{r,h) = h dn^l^^^fA^\/M^. (83) 

Ji du 

r{u), the inverse of u{r), is itself a power series, so that 



dlnr(u) Z'^A" ^ m sr-^ / Cr+i 1 m . 

dli ^^\~h) 1?+^ ^ u ~^ ~h ^ \~h ) IF^ ^ u^ ~h'^^^'' 

s = r = 

(84) 

We have split out the main part 1/u from the correction 0{m/h). C's are 
numbers 0(771") constructed with the set {Nk}: 

Co = 1, Ci = A^i, C2 = iVi + 2N2, Cs = n! + 6N1N2 + 3N3, .... (85) 

Hence 

G(r,/i) = /i^ C.J.(u) = ^7n=G.(r,/i) , (86) 



where 




are elementary functions. Except for constant contributions, their power 
expansions for large u are odd (even) for s even (odd). As implied in Eq. 
(|86[) . h is not expanded in the functions G3. 

With this general formalism we can draw an interesting conclusion 
about enhancement, which corresponds to the limit {u^ ,u^) ^ 1. When 
u^ 1 the functions Js (u) converge to a finite limit of order unity, except 
for Jo{u) u and Ji(ii) Inu; hence, when h is fixed, at higher order no 
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enhanced terms arise in G(r, h) and in the reduced action. Enhancement 
occurs only when h itself is expanded and expressed in in terms of the 
geometric distances of the end points, just as it happens in the case of 
Newtonian hyperbohc motion. 

Using the universal integration variable u, the second-order contribu- 
tion to the hght-time in (|67p has been calculated out with the aid of a 
computer algebra code. We have 



+ —(Ni + 2N2) ( arccos — ± arccos — ) . (88) 
2ho \ rg r^ ' 

With the help of (|70|l and (|77I78I) . this expression reduces to 



^ ° 2r^''B(?-B -'r^cos$^B)(r^ -r^cos^^g) 2ho ' 



The last term in (|65|) requires the derivatives S'i{ho) and S'Q{ho) 



(89) 



hoir^ - cos$^3)(r^ - r^ cos$^ 

S'oiho) ^ Nil ^ ± ^ 

^HIb 

■Ta '\Tb - Ta cos "1>^b ) (r^ - cos <l>^g ) 

The last term in parentheses in (|65|) is therefore 



(90) 



(91) 



\ S'i{h,,) t^lr^^(r^+r,f(\-co^^^^f 

2 Si; {ho) 2rgr^sin2$^^(r^ - cos )(r^ - cos ) ' 

(92) 

Combining this with the first term in parentheses in (|65p . we obtain 



'^'^""^ 2S^{ho) r,r,(l + cos$,J + 2ho ' 
N?r^^ ^ r^^ (8 - 4/3 + 87 - 3e)$, 



r^rg(l -f cos$ ^^) 4r^r^sin$ 
The light-time to second order (|65p is therefore 
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+ mNi In [(r^ +r^+ r^^ )/{r^ + 



+ 



m^I^f^^l^^-—^]; (93) 
r.r„\ 2 sm<I>.„ 1 + cos $ , ' 



+ 



ts-^A^ ^AB + mN^ In [(r^ + y^r| - hl)/{T^ + ^ - hi)] 

(^^l±lEl^ ^ (94) 

r^r^y 2 sin*^^ l+cos<l'^^y' 

in the obtuse and acute cases, respectively. Remarkably, A2 has the same 
expression. This agrees with the result obtained in iW. 

With the same technique, using (|68p . we have computed also the re- 
duced action at the third order. For good measure, here is the result: 



•53 W = 6fe2(,2 „ fe2)3/2 (2^1 + ^N.N, + 3iV3)r^ + 

- ih\Nl + QNiNi + ANs)r^ + Qh^N^Ni + Nz)) + 

- 3/1' (iVi^ + GA^iAfa + 4Af3)r^ + 6/i^(AfiAf2 + A^s)) • (95) 



8 Enhancement 

Enhancement occurs in the obtuse case when and are both much 
larger than 60 = ^0 = 1, so that (|75|) reduces to 

(i- + I =2iV,. (96) 

As hinted in Sec. 2 for the Newtonian case, it is appropriate to formally 
introduce another infinitesimal parameter 60/-R = 1/-R, where R is the 
harmonic mean of the distances (|10|l . When the ratio T^jr^ is 0{R"), 
as we assume, the n-th order harmonic average + is 0(l/i?"). 
The intermediate case bo ~ <C rs, not discussed here, also shows 
enhancement. For instance, it occurs in a nearly parabolic orbit with a 
small perihelion distance pq, as in the case of a solar probe, for which 
even = ARq has been envisaged. The expansion of 

hi= RNi+0{l/R) (97) 

has only odd terms. One should also note that, as can be seen from Fig. 
1, the angle is fixed by the Euclidean experimental setup and should 
be considered independent of m. In the approximation ho — 1 <^ R, 

$^^^n-2/R + Oil/R') 

is slightly less than vr; the law of cosines has been used here. 

We now proceed to discuss enhancement at the second and third order. 
It is convenient to first review the behaviour of the function G(r, h) (|86p 
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in the limit r/h = 0{R) ^ 1. Replacing p with its expression (|38|) and 
expanding, one gets: 



Go(r,h) = h.h{r/h), Giir, h) = Ni [jo{r/h) + Ji(r//i)] , 
G2{r,h) = ^X{r/h) + I^^J^(r/h)+2Ji{r/h)] + ^J2{r/h). 



Now, when 1 



Ji{u) = -1 +ln(2M) + . 
setting u — r/h, 



J2 M = T - - + 
4 M 



Go(r,/i) 


= r — 


Tih 
- + 




24r» 


Gi(r,/i) 


= iVi 




+ ln 


2r 


G2(r,/i) 


= iVi 




gj.3 





+ . 



2/i 6r3 



We need 



Go.hhir,h) = + Go^hhh{r,h) 

^ r ,^ A. / 1 h \ Ni 



GiMn{r,h)^N,{-^ + ^]^^-^, 



ft2 



TT 



iV2 



TV 

'2h? 



h 
3r3 



2N2) 



4/i2 



In the expression (|67|l of A2 the last term is constructed with G2 (r, /i) 
and is not enhanced. The second term comes from Gi^h{r,ho) = — A*'i 
and, when summed over the end points, contributes to the light-time with 
-2Nihi = -2NIR. Lastly, the first term gives hl/R = N^R. Therefore 
the enhanced part of the second-order contribution to the light-time is 

Aacnh = -iVi' 7? + 0(i?°), (98) 

in agreement with ((Ojl. The second-order terms in the ODP are just the 
enhanced ones. 

In a similar way, we get the enhanced third-order terms. For this 
we need the enhanced part of /12 , to be extracted from (|64l) ; its terms are 
constructed, respectively, with Go,hh, Go,hhh, Gi^hh and G2,h. Using their 
asymptotic expressions above one gets the relation 



24 



The third term prevails, and 

hi = -hiNiR = -Afi'i?^ + 0{R), (99) 

in agreement with the Newtonian case (which corresponds to A^i = 1). 
In the expression (168^ for A3 



TT 4 4) 

-^{N? + 2N2)-3iN! +6NiN2+4Ns){r^+r^) (100) 
the first, third and fourth terms are enhanced, so that finaUy 

A3cnh=iVl'i?'+0(ii). (101) 

Similarly, it turns out that A4cnh oc N^R^ + 0{R^). 

To summarize, the expansion (|12p reads (for the Sun): 



^ = Ai + 2 X 10-^:?° A2 + 4 X 10-^^ f^^^^ 
m &o V 00 

In the obtuse case, when R ^ 60, A3 a descending power of R/bo, be- 
gmnmg with (R/bo)"'^. This is the main enhanced term. It depe nds 
only on the single PPN parameter A''i : one could say, enhancement arises 
due to the long-range component cx 1/r of the index of refraction. Ai, 
typically ~ 10A''i, is the (logarithmically enhanced) term of Q; 



j A3 + ... (102) 



0(1) 



o 



single out the main enhanced contribution. For a given R, the strongest 
possible enhancement occurs when bo = Rq; numeriocally 



— = 10iVi -2 xlO-«Aiv2 +4x10-^2 ("A^ +.... (103) 
m Rq \RqJ 

In a typical configuration, with one station on the Earth, Ra ~ lAU <^ 
r_g, so that R — 2 AU — 400 RQ. The three terms in the expression above 
are about 20,3.2 x 10~^,6.4 x 10"''. For a given accuracy in A''i (or Ai) 
this shows how many terms are needed in the expansion in this extreme 
case. 
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9 Deflection 



In the standard theory, the deflection of the image of a far away source 
is the acute angle S between the asymptotes of the ray. Taking the origin 
of longitudes on the symmetry axis OK through closest approach (Fig. 
4) and using (|53p . the longitude of the outgoing asymptote reads (with 
p = uh) 

n + 5 ^ r°° dr ^ dlnr(u) 1 

= — TT- = h I , = / du — ^ , =. (104) 

2 A rsjp^ - Ji du V^?T7T 

Expanding in powers oivn/h, using l|85[l and separating out the main part, 
't>^=Y.CsIs[^) + , (105) 

s=0 s = l 

where 

du 



Is = 



/i u^+WW^ ~ 1 

are numerical constants and d{log{r{u))/du has been deflned in Eq. (|84 
The total deflection is, explicitly 

S = 2N,- + n [-) + [-) +... 



In the more common isotropic gauge (|82p 

h = bN{b) = b + Nim + N2 — + N3-r^ + 

b b^ 

and so 

5 = 2iVi- + ^iN? + '2N2)^4N! ^my ^ 
b 2 \ b J 

lOiVf + I8N1N2 + I2N3 - 3TTNf - GTriViiVa fm\ 3 



+ - 



(=)%.M,», 



In terms of the PPN coefficients and using the expansion of h, to second 
order we have 

S = ^intlll + ^ (8 - 4/3 + 3. + 87) , (107) 

which agrees with J_3^ ; in general relativity, and using the closest approach 
b, 

/-.r o<,^"^^ (155 -457r)m^ 
5 = 4- + (15^ -32)— + ^ 36r^' (108) 

in agreement to second order with [7]. 

This standard approach, however, is not adequate for astrometric ob- 
servations, which are carried out from a point B at a finite distance . 
In the linear approximation this problem has been solved in |19) . Sec. 
40.3; here we give a general formulation and derive the quadratic term. 
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Figure 4: Deflection measured from a finite distance. A ray from a far away 
source arriving along the direction n* is deflected, and arrives at the observation 
point B from a different direction, with unit vector (eq. (I110|) ) tangent to 
the ray. The deflection angle 5^ is smaller than the asymptotic deflection (5, 
the angle between the asymptotes. The origin of longitudes is taken on the axis 
OK through the closest approach, so that ^tx3 = n /2 + 5 /2. The figure also 
illustrates the meaning of the parameter h. The point J is the intersection of 
the tangent through B with a line through O perpendicular to the asymptote. 
It is easily seen that the distance OJ — rh/ yJr'^N(rY — h"^ so that at great 
distance this distance becomes h, which therefore is just the impact parameter. 
In the Newtonian dynamical model h (|17p is a constant of the motion, with the 
same meaning. 
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Referring to Fig. 4, we need the unit tangent vector n{cf>) in the coun- 
terclockwise direction (increasing (j)) at a generic point (r cos (p, r sin (j>) on 
the ray (for simplicity, on the outgoing branch), expressed in terms of the 
function r{(j)): 

, , (r' cos (ji — r sin 0, r' sin + r cos (/)) 

From IMl) 



r ((p) = r , 

so that at B, the tangent vector is 

= - h?cos(j}g - hsm(j>g, ^ p\ - sin (ji^ + /icos^^) 

= {n^^.n^^). (110) 

With 



cosXb = VPb = l/Wfl. sinxB = \jp\ - h'^ / Pb = - l/'^l. 

it is convenient to introduce the quantity Xb , a function on the ray; in 
the limit m — > 0, since p —> r and /i — > 1, it reduces to (Fig. 1). Then 

"b = (sin(xB - <t>B)^ cos(Xb - </>b))- (111) 

The deflection (5^ is provided by the vector product 

|n* X I = sin 5^ , 

where n* = (sin(5/2), cos((5/2)) is a unit vector along the asymptote of 
the incoming ray. Hence we obtain the exact expression 

Sb=^b-Xb+^- (112) 



Two effects contribute in (|112|) : a local term Xb due to the change in 
the tangent, and a change in the orientation of the outgoing asymptote 
relative to OA. In the case of GAIA and other space astrometric projects 
no images can be obtained near the Sun, so that = 1 AU ~ h and 
there is little enhancement. The data analysis will be truly global, with 
subtle statistics. The expected angular measurements error « 5 x 10~^^ 
is quite below the first-order deflection ~ 4 x 10~* and much larger than 
the second-order term ~ 10~^^; but the fractional difference between 5 
and 5g is not small. With our powerful formalism the derivation of the 
second-order approximation to 5b is straightforward. 

Two limits are noteworthy. When m 0, (j>B tends to and, of 
course, there is no deflection. To recover the standard expression when B 
goes to infinity, note that, using [84} , 

, dlnr(it) 1 m , q{u) 
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therefore to second order 



+ ^^^Z + 2h?^^ • ^^^^^ 



Thus, the deflection reads 



In the hmit ii^ — *■ cxa, this agrees with Eq. (|106p . In terms of and ho, 
this is 



5 !~, TT m mCi\ 



m C2Xb _^ m Ci ( N-iho h^r^ 



) ■ (115) 



For Ug finite, this resuh agrees in first order with |19| . Eq. 40.11. 



10 Conclusion 

With the implementation of optical lasers in deep space, experimental 
gravity will undergo a big leap. The planned mission ASTROD ([23], [22] 
and other papers) will consist in a fieet of three drag-free spacecraft in a 
triangular configuration with semi-major axes of about 1 AU. Although 
no detailed error analysis is available, ranging accuracies of 3 x 10^'^ cm 
or better are expected; with closest approach less than 1 AU, this error is 
comparable with, or smaller than, the second-order gravitational delay. 

Optical interferometry in space will make huge improvements in phase 
measurements possible. The GAME (Gamma Astrometric Measurements 
Experiment) project (see [TS]) consists in a Fizeau interferometer in the 
focal plane of a space telescope to measure the angular separation of stars 
in a narrow field of view near the Sun. The expected accuracy in 7 of 
10"'^ wiU require second-order corrections in the gravitational delay. 

LISA - a planned mission for low frequency gravitational wave detec- 
tion ([14], [12] and many other papers, in particular [20]) - will fiy three 
drag-free spacecraft orbiting at 1 AU at the vertices of an equilateral tri- 
angle with sides L = 5 x 10^^ cm; this fieet will rotate around its centre 
with the period of a year. Three optical interferometers with baseline L 
will operate simultaneously, with an expected sensitivity a^jL k, 10~^^ or 
better. The change in light-time difference between two arms due to the 
solar gravitational delay has the period of six months, in a frequency band 
overwhelmed by the acceleration noise, but it is interesting to evaluate the 
effect. For two vertices A and B, — = 5r ~ L ^ ('"ai^b) ~ 1 ^U. 
In the (now generic) acute case the reduced action (1511) (with the — sign!) 
is of order 

"^AB +'mNi—zJL== ^ 5 X 10"cm-f lO^'cm. 
V'^i - hi 

With the approximation 5r <Sil AU the action reads 
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an expression which can be used directly to obtain all relevant quanti- 
ties. For an estimate, however, it suffices to remark that in the above m- 
expansion each term is smaller than the previous one by 0(m/r^) = 10~*; 
hence for LISA the first-, second- and third-order corrections to the light- 
time are, respectively, of order 10* cm, 10~* cm and 10~^^ cm, corre- 
sponding gravitational wave signals of order 

2 X 10"*, 2 X 10"^^ 2 X 10~^*. 

We did not investigate the consequences of this large, but low-frequency 
signal on the performance of the instrument. 

The puzzle of the ODP expression for the gravitational delay has been 
understood. It must be considered in the framework of an expansion in 
powers of m/h; of all second-order terms so arising, in a close conjunction 
some are enhanced. They can be rigorously singled out with a further 
expansion in diminishing powers of R/bo; those that appear in the ODP 
are just those of order m{m/ bio) (R/bo). With the powerful tool of geo- 
metrical optics, we have provided a procedure to extend the calculation 
to higher order and have obtained the full correct second-order term of 
the delay. 

A methodological reflection is a flt conclusion. The evaluation of the 
gravitational delay, a conceptually simple and straightforward problem, 
faces subtle mathematical difliculties and a great algebraic complexity. 
Our approach is based upon two unusual mathematical levels of descrip- 
tion: light propagation with the eikonal theory, rather than null geodesies, 
and asymptotic power series, an abstract mathematical tool. The latter, 
in which ordinary functions are set aside and an abstract mathematical 
tool is employed, seemingly runs against physical intuition. As shown, 
both are essential to directly attain, and take advantage of, the crucial 
features of the problem: the light-time as the minimum of Fermat's ac- 
tion, and a safe and automatic procedure to select and estimate different 
terms. This is another example of the tenet that every physical problem 
has an appropriate, often not intuitive, level of mathematical description, 
and severe penalties are in store for its neglect. 

Appendix 

The radial gauge freedom and the difference between closest approach and 
bo can cause some confusion. For example, the textbook [5D] presents (eq. 
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(8.7.4)) the light-time between closest approach and a generic point; it is 
expressed in Schwarzschild's gauge r and reads 



quite different than ([2]). In a real case two such terms are needed, one 
for each branch. But, contrary to what stated in the textbook, the sum 
of the two square roots (first term) is not the distance AB. The isotropic 
gauge and the distance 60, not the closest approach, should be used. First, 
setting (|37|) f = r + 7m, the formula reads, to 0(m), 



Both formulas are useless, however, because the closest approach 6 = 
bo + mfoi = 1 + mbi is not known beforehand. The ray must be anchored 
to two known points and, with 61, is determined by the unknown 7 with 
([801). Since 





\/r2 — 62 _ ^/r^ — 1 



6: 



m 




and the standard formula is recovered. 
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List 



of symbols 



A event or point where the photon starts 

A{r) metric coefficient 

B event or point where the photon is detected 

B{r) metric coefficient 

b closest approach in isotropic variable 

C(r) metric coefficient 

h closest approach in Moyers's variable 

60 Euclidian approximation of the same 

£ Euclidian arc length 

m gravitational radius 

N{r) — -^pjy index of refraction 

ODP Orbit Determination Program 

Pq perihelion distance 

— ^ A^B harmonic mean of the distances 

r isotropic radial coordinate 

_Rq radius of the Sun 

r(/) photon trajectory 

P0 perihelion distance 

S Fermat's action 

Q(x^) eikonal function 

t time in the rest frame of the mass 

tj^ starting time of photon 

arrival time of photon 

7 relativistic PEN coefficient 

gravitational delay 

As expansion coefficients of delay (|12p 

A undefined parameter along the light path 

p = rN{r) Moyer's radial coordinate 

(j> longitude 

Phi longitude 
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